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Abstract 

Aedes aegypti is known as the responsible vector in transmitting 
dengue flavivirus. Unavailability of medication to cure the transmis¬ 
sion in human blood becomes a global health issue in recent decades. 
World epidemiologists are encouraged to focus on inves- tigation to¬ 
ward an effective and inexpensive way to prevent dengue transmission, 
i.e. mosquito control. In this paper, we present a model depicting the 
dynamics of mosquito population based on indoor?outdoor life cy¬ 
cle classification. The basic mosquito offspring number is obtained 
and analysis of equilibria is shown. We bring along a discussion on 
application of optimal control model which engineers two simultane¬ 
ous schemes. The first scheme is done by disseminating chemical like 
Temephos in spots where eggs and larvae develop, meanwhile the sec¬ 
ond scheme is done by conducting fumigation through areas where 
adult mosquitoes prevalently nest, indoor as well as outdoor. A ver¬ 
sion of gradient-based method is presented to set down a workflow 
in minimizing the objective functional with respect to control vari¬ 
able. Numerical results from analysis of the basic mosquito offspring 
number with constant control and from optimal control suggest that 
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one has to enhance the usage of fumigation rather than Temephos. It 
is also suggested that applying both control schemes simultaneously 
gives the most significant reduction to the population. 


Keywords: mosquito population dynamics; mosquito reduction man¬ 
agement; the basic mosquito offspring number; optimal control. 


MSC: 49J15; 49M05; 92D30. 

1 Introduction 

It is known that Aedes aegypti mosquitoes be the primary transmitters of 
dengue fever in the world [2], This species were predominantly tropical 
species and confined to coastal areas. Now they are widespread inland and 
cause deadly morbidities by means of dengue fever in mostly Southeast Asia, 
Africa, USA, Australia, Brazil, Argentina, Caribbean Islands, China, India, 
Japan, Portugal, etc mm- The areas where Aedes aegypti mosquitoes most 
probably prevalent in which have tremendous threat of dengue outbreaks are 
densely urbanized areas M- Several studies confirmed that eradicating 
mosquito population had constricted the area of dengue endemicity through¬ 
out the globe, beside it helped to bring down the cost as regards a high 
number of surveillance. Even vector control has been considerably inexpen¬ 
sive, a government tendency to take out of concern the funding to attempt 
it continuously leads to a necessity of a more efficient and effective control 
management. In order the mosquito population to be precisely managed, one 
needs to model the dynamics of such population and further to investigate 
effectivity of a control scheme which is acted to it. 

Aedes aegypti belongs to species which inhabit in domestic water contain¬ 
ers: bath vessels, flowerpots, drums, tins, unused tyres, untreated swimming 
pools, or even in curved broads where it is possible for water to last on for 
long time HI. Both male and female adult mosquitoes feed on nectar. Only 
females require additional blood sources to obtain nutrients before producing 
eggs and during eggs maturation. Note that a single female can lay down at 
one time about 100-200 immersing eggs in water [5j. Potential eggs can be 
produced by Aedes aegypti females up to five batches during their lifetime. 
An egg needs 2-5 days to maturate its living embryo, depending on water 
temperature mm- In an advanced growth, each egg turns to a larva which 
posteriorly withstands for 5-10 days depending on water temperature and 
air humidity HU- This living larva usually eats algae and microorganisms in 
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water surface, making each individual competes with its owns for logistics. 
A discussion on competitions amongst larvae can become important after 
the fact that in a joint container, Aedes albopictus larvae outcompete Aedes 
aegypti larvae, thus the winners develop at a faster rate [1]. In the prescribed 
range of lifetime, each larva undergoes four times skin exfoliation and turns 
to a pupa at the very end of the processes. An idle pupa needs to wait for 
1-5 days before it metamorphoses into an adult [16]. In addition, the living 
adult can generally survive for 10 days, or in some extreme cases for 2-4 
weeks [23] . 

Tabachnick et al. m found three polytypic origins where Aedes aegypti 
breeds: domestic (urban housing including its narrowing environment), syl¬ 
van (rural areas, some cases in tree holes and leaf axils), and peridomestic 
(artificial plantation areas such as coconut groves and farms). In general, 
separation of the origins as indoor and outdoor had also been highlighted. 
In [12] , the authors notified that the number of Aedes aegypti differed based 
on indoor-outdoor classification and heterogeneity of containers. A brief cor¬ 
responding result showed that Aedes aegypti species constituted as the most 
abundant in indoor containers compared to the other tested species. Some 
reference also mentioned that in indoor containers, competition amongst lar¬ 
vae had not always been the case since a particular tendency made Aedes 
aegypti grew more than Aedes albopictus |5|. 

By a basic idea of incorporating control measures to a mathematical 
model of mosquito population dynamics, several control schemes were tested 
toward fighting the spread of Aedes aegypti. We highlight, for instances, 
utilization of ultra-low-volume (ULV) insecticide [15] also combined with 
Temephos dissemination [ 22] and the sterile insect technique (SIT) [3 HE] 
(2D HD] as a Genetic-based Vector Control (GVC). However, in this paper we 
do not take the epidemiology of mosquito population or infection exposures 
amongst individuals into account. Therefore, a susceptible-infected segre¬ 
gation is no longer in use. We accentuate a model of mosquito population 
dynamics based on indoor-outdoor life cycle classification and introduction 
of control intervention as well as investigation to the best strategy for re¬ 
ducing the population with cost as cheap as possible. As control devices, we 
highlight utilization of Temephos (mainly used to kill larvae) and fumigation 
(used to kill adults by blocking respiration process) after the fact that they 
have already been well-known schemes in an integrated mosquito eradication 
programme. 

We organize the rest of the paper as the following. In Section 2, we set 
up a mathematical model for depicting the dynamics of multi-age-segregated 
mosquito population based on indoor-outdoor life cycle classification. We 
add two control measures to the system in order to check how the popu- 
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lation trajectories response to such control. In Section 3, we examine the 
biological meaningfulness of the model based on positivity of the system’s 
solution and existence and stability of equilibria. So far we use plausible 
constant value for the control. In Section 4, we discuss the optimal control 
model which generates a time-variant control solving the proposed optimiza¬ 
tion problem. In seeking an optimal solution, we first use the so-called in¬ 
direct method to generate state-adjoint-gradient systems of equations and 
then utilize a gradient-based method to solve the generated systems in an 
algorithmic workflow. As the rest of the paper, in Section 5, we do numerical 
tests for bringing forward some visualizations of the model. 

2 The mathematical model 

To capture the dynamics of mosquito population, let us first classify the 
population based on age-segregation: egg E. larva L and adult A. We out- 
concern the dynamics of pupa by our assumption that there will be no deaths 
and imbalance in- and outflow occurring amongst pupae, since then the pop¬ 
ulation of pupa remains constant. Based on indoor-outdoor classification, 
the eggs and larvae differ from indoor Ei,Li and outdoor E 2 ,L 2 . No more 
segregation for the adults since every individual can fly to wherever it preva- 
lents, indoor as well as outdoor. We introduce the control measures u x and 
u 2 denoting the impact rate of utilization of Temephos dissemination and 
fumigation respectively. Let [0, T] be the range of observation time t and 
x = (Ei, E 2 , Li, L 2 , A) t e L 2 ([0, T ], IR 5 ) be the state variable. Our model is 
preliminarily exhibited by the following equation 

x(f) = f(x(t), u(t)), x(0) = x 0 y 0. (1) 

We assume that adult mosquitoes will select indoor breeding sites with 
the probability p and, therefore, select outdoor breeding sites with the prob¬ 
ability 1 — p. Alongside with introduction of the tendency-based probability, 
we denote p as the corresponding rate of the adults to lay eggs. It is as¬ 
sumed that natural deaths can occur in all classes, thus we denote , 5 } as 
the corresponding rates. In average periods, both indoor and outdoor eggs 
metamorphose into corresponding larvae with the transition rates aq x , 2 }- The 
same situations hold for respectively indoor and outdoor larvae, that they 
metamorphose into adults with the transition rates @{ 1 , 2 }- As an underlying 
discussion, we introduce the two logistic coefficients cq Xi2 } to accommodate 
phenomena of competition amongst larvae. Assume that all in- or outdoor 
water containers are treated as homogeneous such that — proportionates to 
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some constant M. The fact that suitable choice for M > 1 leads to the sit¬ 
uation where outdoor breading sites have M-times larger carrying capacity 
than indoor. After all, the model (JT)) is unfolded as 


ppA - («i + r/ijEi - quiE 1 

(2a) 

/i(l -p) A - (a 2 + 77 : 2 )E 2 

(2b) 

opEx — oyL 2 — (/?! + r/ 3 )L| — uxLi 

(2c) 

CV2E2 — O2L2 — (/?2 + ? i 4 )Lo 

(2d) 

ft lL X + ft2^-‘2 — Vo A — u 2 A. 

(2e) 


Two control measures are added into the model. Temephos is practically dis¬ 
seminated into indoor water containers to kill larvae and eggs. The weighting 
factor q G [0, accounts for the fact that the dissemination has less impact 
on eggs. Meanwhile fumigation is directly targeted to adult mosquito popu¬ 
lation. 

It is assumed that all the control u belongs to a set of admissible control 
U C L 2 ([0, T], 1R+). Define 

V := {x : x(t) = f(x(t), u(t)), t G [0, T], x(0) =x 0 ^O,u6W} (3) 


as the feasible set of solutions. It is intuitively believed that the higher u 
plugged into the model, the lesser the number of all classes. On the other 
hand, the lesser u, the more the population tends to explosion. The more u 
means that a policy maker needs to spend more fund, otherwise there will 
be no significant reduction to the mosquito population size. Given positive 
trade-off constants w x ,{i,-, 5 } and w u ,{i, 2 }, the following objective functional 
accommodates the necessity of balancing situation between significant pop¬ 
ulation reduction and limitation of fund: 

2 Tj 0 

Now our optimization problem is stated as 



5 2 

w x,;Xi (0 + J2 d u ? (*) dt • ( 4 ) 


(OC): find (x, u ) eT> xU such that J(u) —» min . (5) 


3 Model analysis 

We need to ensure that our model is biologically meaningful. The following 
theorem gives a primary meaningfulness of the model: whenever the initial 
condition is positive, then evolution of the solution points in forward time 
stays positive. 
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Theorem 1 Consider the model (0Q) where uEU. If xq >z 0 then x(t) y 0 
for all t > 0, or particularly on the observation range (0, T], 


Proof Following the steps in [22], we let n be (5 x 5)-matrix representing 
a collection of all normal vectors (by rows) to the boundary of nonnegative 
orthant cdR/)_. Thus we have n = —15 where 1 5 denotes the identity matrix. 
To ensure that the trajectory of solution does not walk out of the nonnegative 
orthant, one only needs to check the trajectory of points in the boundary. 
Notice that at i-th boundary, djlR/L, 


[ nf ( x i u )]JxediiR5 ,uew < °- 

This means the direction of evolution of such points in a boundary is in 
counter-direction or at least perpendicular to the corresponding normal vec¬ 
tor. Thus it follows that the solution must not leave IR+ for all t > 0. □ 

Let us first consider the autonomous system of nonlinear differential equa¬ 
tions ([I]) with constant control u G IR+. For abbreviate purpose, let s 4 = pp, 
s 2 = p(l - p), s 3 = l, s 4 = a 2l s 5 = Pi, Sg = p 2 and di = aq + rp + qui, 
d 2 = a 2 + rj 2 , d 3 = A + rj 3 + u 4 , d A = p 2 + rp, d 5 = p 5 + u 2 . In order to 
obtain equilibria of the system, we need to solve f(x, u) = 0. It follows from 
this process, one holds 


_ s 5 L] + sgL.j 

G?5 

g* = Sj(s 5 Lj + s 6 Lj) 
did§ 

* _ s 2 (s 5 Lj + s 6 Lj) 
2 d 2 d 5 

where (Lj.L^) follows the leading equations 

Lj T GqL 4 T 6 4 L 2 = 0 
L.j + a 2 L 2 + & 2 Li = 0. 


( 6 a) 

( 6 b) 

( 6 c) 

( 6 d) 

( 6 e) 


All the constants that belong to the last equations are given by 

d 3 S1S3S5 S1S3S6 d 4 s 2 s 4 Sg 

°i —- ~r~j —j h — --T-. —, 0-2 —- -j —,—, b 2 — - 

(J 1 did 5 ai did^ax a 2 d 2 d 5 a 2 


■S 2 S 4 S 5 

d 2 d$o 2 


It follows from (l 6 d|h (l 6 cj) . three equilibria can simply be obtained 


Qi — 

Q 2 

Q; — 


( 0 , 0 , 0 , 0 , 0 ) 

SlSgOi s 2 S5a 4 


d\d^ 

SlSgCl 2 

d\d§ 


d 2 d§ 

s 2 s 6 a 2 

d 2 d§ 


, n 4 ,0, 
, 0, —a 2 , 


S 5 U 1 

A 

Sq 0 2 

dr, 
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Definition of Q 2 and Q 3 only makes sense in biological point of view if only 
ai,d 2 are negative. The definition simply means that after very long time, 
the class of either in- or outdoor larvae tends to extinction while the other 
classes stay alive. In biological context, such situation can hardly happen. 
In the next writing, we demonstrate that the choice of ai,a 2 to be positive 
leads to a more interesting discussion. First, we introduce a measure whose 
cubic value is given by 


n «r - 


3_ S4.S3.S5 ■S2S4S6 


+ 


d$did 3 d 3 d- 2 d± 

The following identity is found after some algebraic computations 

6162 — O1O2 


d 3 d. 


3«4 


cricr 2 = 1Z{vlY - 1. 


(7) 


( 8 ) 


Using unity as a threshold value, the following theorem justifies the stability 
of zero equilibrium Q 1 by considering the nominal of 77(u) relative to the 
threshold. Meanwhile, the next theorem shows the existence of a positive 
nontrivial equilibrium (also well-known as coexistence equilibrium). 

Theorem 2 The zero equilibrium Q l is locally asymptotically stable iflZ(u) < 
1 and is unstable iflZ(u) > 1. 


Proof Checking the stability of the zero equilibrium in local view, it is 
similar to see behavior of solution of the linearized system (jTJ) in the neigh- 


_,u)(x 

-Qi 

). In 

this case, 

L . Unfolding 

this 

A Qi, 

u), 

—di 

0 

0 

0 

Si 

0 

-d, 2 

0 

0 

S 2 

S3 

0 

CO 

1 

0 

0 

0 

S 4 

0 

—d 4 

0 

0 

0 

S5 

s 6 

~d. 


<7(Qi,u) = 


Q x is locally asymptotically stable if and only if all eigenvalues of J(Q l5 u) 
have negative real part. To see this, a simple cofactorization method com¬ 
putes determinant of J(Q 1; u) — A/5 as the sum of all its cofactors, given by 
negative of 

(di+A) (d 2 +A) (c/j+A) (d 4 +A) (ds+A)—s 2 S4S6(di+A) (c/H-A)—siS3S5(d 2 +A) (</ 4 +A). 
Clearly coeff (A^ 3 , 4 ’ 5 1 ) are real positive, meanwhile 
coeff(A°) = d\d 2 d 3 did 3 — sis^s^d^ — s^s^s^didj, 

coeff (A 1 ) — d\d 2 d 3 d 4 T ■ ■ ■ T d 2 d 3 d^d^ — 545355(^2 T d 4 ) — 'S 2 S4Sg(di T dY) 
coeff(A 2 ) = d\d 2 d 3 H-h d 3 d 4 d 5 - S4S 3 s 5 - s 2 s 4 s 6 . 
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If 7£(u) < 1, or equivalently ^^5 _[_<- ^ one can see that coeff (A^ 0,1,2,3,4,5 ^) 

are real positive. Since lim A ^. ±00 — det(J(Q 1 , u)—XI 5 ) = ± 00 , — det(J(Q 1 ,u) — AJ 5 )Ia=o > 

0 and 


ddet( J(Q 1 , u) 

dA - 


a/ 5 ) 


> 0 , 

A >0 


d 2 det(J(Q 1 , u) - AJ 5 ) 
d\ 2 


> 0 , 

A >0 


then det( J(Q 1; u) — A/ 5 ) can not have zero that has positive real part. More¬ 
over, all zeros of det(J(Q 1 ; u) — A/ 5 ) have negative real part. If P(u) > 1, 
then coeff (A 0 ) < 0. Then for however the values of coeff (A 1 ) and coeff (A 2 ), 
it follows that there must be at least one real positive zero. This means Q 1 is 
unstable. If P(u) = 1 , then a similar step can prove that det( J(Q 1 , u) — AJ 5 ) 
has one zero with the value 0 and the other zeros remain with negative real 
part. In this case, Q 1 is considerably still stable but not asymptotical. □ 


Theorem 3 If a\ and 02 are positive, a unique coexistence equilibrium Q 4 
exists if IZ(u) > 1 and does not if 7Z(u) < 1. If a\ and 02 are negative, 
then there always exists a coexistence equilibrium for any value ofTZ(u). If 
else, then it always holds IZ(u) > 1, and therefore, there always exists a 
coexistence equilibrium. 


Proof One simple substitution in (I6d|h(1(3c]l makes the equilibrium state L 2 
to follow 


Li{ + 2 oiL^ + (a 2 — 6io 2 )lT T (6 2 6 2 — a.i&ifl2)Li — 0 . 

As Li = 0 results to equilibria that has been discussed before, now let us 
consider the remaining cubic 

-P(Lu) := Lj* + 2 oiLj + (a 2 — 6 io 2 )Li + (6 2 6 2 — 0161O2). 

Since coefficient of L) is 1 > 0, then liniL^ioo P = ± 00 . We prove this 
theorem by parity. 

Case 1 Oi and 02 are positive. In this case, three possibilities are carried 
out: 7?.(u) > 1, 7?.(u) < 1 and 7?.(u) = 1. If 7\l(u) > 1, then by identity (JS]) it 
holds that 6162 > 0102 - Thus we have 0 > b\b 2 — 016102 = P(0) since 61 < 0. 
It follows from the value of P( 0) and liniLi-i-oo P = 00 that P must have at 
least one positive root. Now we have to check that L 2 = Li(L ) 1 b ) l ~ ai >0. To 
see this, let Li = p 2 be the corresponding root, where p G 1R\{0}. Extracting 
(Li — p 2 ) out of the cubic, we get L 2 + (2ai + p 2 )Li + J\f where 


J\f :— (o 2 — 61 O 2 ) + p'(2ai + p 2 ). 








Since J\f = ^ lbia ' 2 ' 1 > 0 then together with the claim of uniqueness of 

positive root L | = p 2 , we need (2a! + p 2 ) > 0 to ensure that the remaining 
quadratic function does not have real positive root. Since a 2 > 0, then 
(p 2 + a\) >0. Now we have proved that L 2 is positive, and by ([6]), the other 
states are also positive. If 72.(u) < 1, then P(0) > 0. Since liniL^-oo P — 
—oo, this means there exists at least one negative root Li = — e 2 of the cubic, 
where e E 1R\{0}. Now extract (Li + e 2 ) out of the cubic and we get the 
remaining quadratic function L 2 + (2ai — e 2 )L! + A4 where 

M := (a 2 — bia 2 ) — e 2 (2ai — e 2 ). 

It turns out that Ai = b i b2 ~^ bia2 — LA > o. Thus further analysis confirms 
that a positive root exists whenever (2ai — £ 2 ) is strictly less than zero and 
the discriminant A = (2ai — e 2 ) 2 — 4A4 > 0. One can see that A = (2ai — 
e 2 ) 2 — 4(a 2 — & 2 a 2 ) + 4e 2 (2a 2 — e 2 ) = (2a! + 3e 2 )(2a! — e 2 ) — 4 (a 2 — Aa 2 ) < 0 
since ai, a 2 > 0 and b\ < 0. This contradicts the necessity that A > 0. Then 
there does not exist positive root of the cubic. At last, if 72. (u) = 1, then 
b\b 2 — aib\a 2 = 0 and we get the remaining quadratic L 2 + 2aiLi + (a 2 — b\a 2 ). 
Since ai, a 2 >0 and b\ < 0, it is clear that the quadratic does not have any 
real positive root. 

Case 2 a\ and a 2 are negative. If 7£(u) > 1, then it is easily seen that one 
positive root exists. Now we have to check that L 2 = Ll( ' L . 1 b ( hai ^ > 0. To see 
this, as it is previously done together with the claim of uniqueness of positive 
root Li = p 2 , we need (2ai + p 2 ) > 0 to ensure that the remaining quadratic 
function does not have real positive root. Then (p 2 + ai) > (2ai + p 2 ) > 0 
since ai < 0. Now we have proved that L 2 is positive, and therefore, the 
other states are also positive. If 7£(u) = 1, we have b\b 2 — aib\a 2 = 0 
and the remaining quadratic function is L 2 + 2a 2 Li + (a 2 — b\a 2 ). Since 
A = 4a 2 — 4(a 2 — b\a 2 ) = 4bia 2 > 0 and 2a 2 < 0, then a positive root exists 
whose uniqueness is confirmed by the value of (a 2 — b\a 2 ). Thus we need to 
proof that L 2 = Ll( ~ Ll 6 ^ ai) > 0. Let L 2 = p 2 be such root, then the other root 

must be Li = —(2ai+p 2 ) where (2ai+p 2 ) = — ^ ai ^ a2) . However, the positive 

root is unique if (a 2 — b\a 2 ) < 0, then p 2 + a 2 > p 2 + 2ai = _> o. 
This conhrms that L 2 > 0. If 72.(0) < 1, then as before, we have at least one 
negative root L x = — e 2 and the remaining quadratic L 2 + (2a 2 — e 2 )Li +JH. 
Observe that M. > 0 and (2 a x — e 2 ) <0 since a 2 < 0. Some factorization 
confirms that there exist two remaining positive roots l[ 1,2 ^ = {p 2 , £ 2 } where 
£ 2 = — (2a 2 — e 2 + p 2 ) and p 2 ^ 2 = A4. Then we have to check the result 
of L) ’ + a since this conhrms positivity of L 2 , or furthermore, existence of 
positive equilibrium. Since (p 2 + ai) + (£ 2 + a 2 ) = — (2a 2 — £ 2 ) + 2a 2 = e 2 > 0, 
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this ensures that at least one of the roots is greater than —ai, since then at 
least one root makes L 2 positive. 

Case 3 Else. Since it always holds > 01 ^ 2 , then this implies 7£(u) > 1 
and P(0) = b\b 2 — ai&ia 2 < 0. Therefore, there is always (not necessarily 
unique) coexistence equilibrium. This completes the proof. □ 

All classes have responsibility to generate offsprings in sense that new¬ 
borns are impossible whenever one class in the system remains zero for all 
time. In line with this sense, we define the mosquito-free equilibrium (MFE) 
as a static condition when all classes can not reproduce for the next offsprings. 
This means MFE is equivalent to the zero equilibrium Q x . Following the next 
generation method [8], we define two matrices F(u),V(u) where F(u) is a 
matrix with zeros in the main diagonal and V (u) is a positive diagonal ma¬ 
trix such that J(Q 1 ,u) = F(u) — V(u). Then we define the so-called next 
generation matrix G(u) = F(u)C(u)~ 1 which is nothing but 


G(u) 


( 0 0 0 0 t \ 

0 0 0 0 g 

g 0 0 0 0 

0 | 0 0 0 

V 0 0 f f 0 / 


(9) 


(i, j)-th element of this matrix represents the average number of new indi¬ 
viduals in compartment i produced by a single individual from compartment 
j during compartment j's average individual lifetime period. It is easy to see 
that the spectral radius of G(u) is given by 


max{|p| : u G R+, det(G(u) — plf) = 0} = TZ( u). (10) 


This TZ( u) is often called as the basic mosquito offspring number. Formula¬ 
tion in (U) provides dependency of the basic mosquito offspring number on 
several parameters and the control u. In Section 5, we check behaviour of the 
basic mosquito offspring number with respect to the control and to several 
unknown empirical parameters. 

Theorem 4 The existing coexistence equilibrium Q 4 is globally asymptoti¬ 
cally stable in the nonnegative orthant 1R^_ if IZ(u) > 1. 


Remark We omit writing the detail for the proof of the theorem. Further¬ 
more, it requires the use of the so-called Krasovskii -LaSalle principle. As an¬ 
other reference, one can further see from numerical results that the solution 
of the model is monotonic and tends asymptotically to Q 4 , by supplemented 
condition that the set of parameters is chosen such that 7£(u) > 1. 
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4 Optimal control problem 


Recall our optimization problem: 

find (x, u) G V x U such that J(u) —» min . (OC) 

Denote by U a set of admissible controls where it is assumed to be compact. 
Let g be a functional such that J(u) = f 9 g(x, u) d t. Let u be an optimal 
control that solves (IQCD. Consider a small variation around the optimal 
control 

u e (t) = u(t) + e/c(f), |e| 1. (11) 

Plugging this variation into (Jl]) together with positive initial condition, we 
have the resulting perturbed state x e , where it holds 

x e — x = f(x e , iT) — f(x, u). (12) 

Note that f G C' 1 (1R 5 x IR 2 ,IR 5 ), thus we have the following Taylor approach 

-^(x e - x) = f x (x, u)(x e - x) + ef u (x, u)« + 0(e 2 ). (13) 

Working at 0(e 2 ), we let cp be the solution of the differential equation 

(p = f x (x,u)^ + f u (x,u)«, p(0) = 0. (14) 


Thus some algebraic computations show that our perturbed state is given by 


x e = x + eip + 0(e 2 ). 


(15) 


Further it can be proved that the solution of (1141) exists and therefore the 
term exists. 

Let z G L 2 ([0, T], 1R 5 ) be some dual variable. It follows that -^(z ■ 
x e ) d t = z(T) ■ x e (T) — z(0) • x 0 . Append this to the objective functional and 
we get the following expression 


J(W) 



g(x e ,iT) + 


d 
d t 


z ■ x 


dt + z(0) • x 0 




(16) 


Since u is a local minimizer of J, then derivative of J over e where e = 0 exists 
and equals to zero, in other words 93 - = 0. Steering (TT6l) together with 

e=0 

simplifying the result by factorization, we get 

dx. e 


0= / (g x + f x z + z) 


de 


/ (g u 


f„z) • ndt ■ 


z(T) 


dx e 


e=0 


de 
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in 

e=0 

(17) 












Note from (1141) - (1 1 5[) that variation of k makes varying. Zeroing the right 
hand side of (fT71) together with taking k = g u + f^z, we get the following 
theorem. 

Theorem 5 Consider the optimization problem (IQCD. Let u Eld be a local 
minimizer for J and x E T> be the resulting state. Then there exists a dual 
variable z E L 2 {\ 0, T], IR 5 ) such that the tuple ( x, u, z) satisfies the following 
system 


dTL 

x = —— with cc(0) = Xq >z 0, 
oz 


z = 


m 

dx ’ 




z(T) = 0 (18) 


for all t E [0, T\. The functional TL(x, u, z) := g{x, u) + (z, fix, u )) is Hamil¬ 
tonian functional, meanwhile all equations in (ITS]) are respectively the state, 
adjoint, gradient equations and transversality condition. 

The adjoint (with transversality condition) and gradient equations can now 
be unfolded as 


and 


Zi = 

—1 w x ,iEi + 

(a 4 + gu 4 + 77 4 )zi - CK 1Z3, 

zi(T)=0 

(19a) 

z 2 = 

-W x ,2 Eo + 

(ck 2 + 1 / 2 ) z 2 — a 2 z 4 , 

z 2 (T) = 0 

(19b) 

z 3 = 

— Wx,3Ll + 

(2 ctiL 4 + fix + u 4 + %)z3 — /3 4 z 5 , 

z 3 (T) = 0 

(19c) 

z 4 = 

—cj Xj4 L 2 + 

(2a 2 L 2 + /3 2 + ^ 4 )z 4 — /3 2 z 5 , 

z 4 (T) = 0 

(19d) 

z 5 = 

— Wx, 5A + 1 

(u 2 + ^5)2.5 - /h» z i - h(! - P) z 2, 

z 5 CO = 0 

(19e) 



w Ul iU 4 - gE 4 zi - L3Z3 = 0 


(20a) 



w U)2 u 2 - Az 5 = 0. 


(20b) 


It is essential that U is bounded as this supports meaningfulness in practical 
implementation of the control. For the sake of simplicity, we consider U = 
L 2 ([0,T], [v!, w x] x [v 2 , w 2 ]) and the according projection 


sat(u) = max (v, min (w, u)) 


( 21 ) 


mapping any control u into U. 

To solve the optimal control problem (IQCD for bounded control u E IA, 
we first consider the unconstrained convex optimization problem 

find (x, u ) E L 2 x L 2 such that J(u) —> min 

and use the projection (ITT!) to map the control into a convex set IA. Since IA is 
a constriction of L 2 , the projection leads us to the optimal solution of (IQCI) . 
After all, Algorithm 1 illustrates our scheme to solve (IQCD . 
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Algorithm 1 Gradient-based method for solving (IOCI) . 

Require: k = 0, an initial guess for the control u fc G U , an error tolerance 
e > 0 and an initial step-length A > 0. 

1: Compute x(j 4— x fc (-;u fc ) and z(j 4— z fc (-; x A '(-; u fc )) consecutively from 
the state (with forward scheme) and adjoint equation (with backward 
scheme). 

2: Compute the objective functional J(u k ). 

3: Compute u fc (x(j,z(j) from the gradient equation. 

4: Update u fc+1 (A) 4— u k + Au fc and u fc+1 4— sat(u fc+1 ). Compute x(j +1 and 

yk~\~l 

Z ll 

5: Compute J(u fc+1 ) and set A J 4— J(u fc+1 ) — J(u fc ). 

6 : if |AJ| < e then 

7: Stop. Set (x, c, J) (x(j +1 , u fc+1 , J(u fc+1 )). 

8 : end if 

9: while A J > 0 do 

10: Update new A 4— argmin se [ 0; A] 0(s) := J(u fc+1 (s)) where 0 is a 

quadratic representative of J with respect to the step-length s. Note 
that the solution exists since 0(0), 0 7 (O) and 0(A) can be computed di¬ 
rectly. Compute new u fc+1 (A) u k + Au fc and set u fc+1 sat(u fc+1 ). 
Then compute new x(j +1 and z(j +1 . 

11 : Compute J(u fc+1 ) and set A J -s— J(u fc+1 ) — J(u fc ). 

12 : if |AJ| < e then 

13: Stop. Set (x, u, J) 4— (x^ +1 ,u fc+1 , J(u fc+1 )). 

14: end if 

15: end while 

16: Set k 4— k + 1. 

17: Go to Step 3. 

18: return The tuple (x, u, J). 


5 Numerical results 

Several discussions in this section aims at determining the impact of control 
intervention to the abundance of mosquito population in visual statements. 
We accentuate in this paper that significant reduction of the population 
evolution in forward time be a carry-over from only a cheap optimal control. 
Realistically, optimal control still suffers from practical drawbacks that one 
can not easily implement it in real situation. The reason is slightly supported 
by undisciplined executors in the field to implement the programme as it 
looks as in optimal control result. On the other hand, a constant control 
benefits from its easy-to-implement scheme as the executors only need to 
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deal with fixed-fund allotment problem with flat daily distribution. Another 
benefit from using constant control is that nominal of the basic mosquito 
offspring number can explicitly be computed. This computation educes some 
important statements regarding endemicity of the observed area. It has been 
preliminarily known that whenever the basic mosquito offspring number is 
less than unity, then two facts arise: the zero equilibrium is stable but the 
coexistence equilibrium does not exist. If this number is greater than unity, 
then the zero equilibrium is unstable and the coexistence equilibrium exists 
and is stable. This roughly means that the endemic situation will emerge 
and stay uninterruptible. Otherwise, constant control result always emerges 
with a higher cost compared to optimal control. 

In order to do numerical tests, all the parameters involved in the model 
have to be represented in numerics. Table [T] gives an estimate for all param¬ 
eters. 


Parana. 

M 

P 

h 

hi 

V2 

h3 

h4 

h5 

<1 

Unit 

- 

- 

day -1 

day -1 

day -1 

day -1 

day -1 

day -1 

- 

Est. val. 

2 

0.4 

3.1 

0.02 

0.01 

0.002 

0.01 

0.4 

0.04 


Parana. 

oy 

ay 

CTl 

o-2 

ft 

ft 

Unit 


“d^ 

inch -1 x day -1 

inch -1 x day -1 



Est. val. 

0.15 

0.13 

0.004 

<71 

M 

0.08 

0.05 


Parana. 

T 

W X ,{1,..., 4 } 

^x,5 

W U ,{1,2} 

V 

w 

Unit 

day 

- 

- 

- 

day" 1 

day -1 

Est. val. 

150 

1 

2 

4 x 10 4 

0 

[1,1] 


Table 1: Set of all parameters used in the model. 


5.1 Constant control 

To some reasons, the basic mosquito offspring number 7ft u) reflects the hier¬ 
archy of endemicity of dengue. The higher 7ftu), the more rapid mosquito’s 
growth, and therefore the higher the number of dengue incidences. In this 
case, we demonstrate primary impacts of using constant control to 7ftu). 
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Alongside with this illustration, we show a dependence of H(0) to some un¬ 
known parameters that may be hard to find, or at least hard to estimate. 
Such parameters can be like p (the prevalence probability of mosquito adult 
to breed indoor) and p (the percapita birth rate of eggs). The following 
figures give the illustrations. 




Figure 1: Contour of the basic Figure 2: Contour of the basic 
mosquito offspring number 7l(u) in mosquito offspring number 71(0) in 
(ui, Ui)-plane. (p,//)-plane. 

From Fig. |TJ some conclusions are highlighted. It is clear from the figure 
that u 2 gives a more significant decrease to nominal of the basic mosquito 
offspring number compared to iq. In constant control case, each value of 
control taken in the set of real positive numbers can roughly represent a 
negative return. A negative return ui of 1 simply means that, 100% the 
number of individuals as in the previous step (in day) has to be killed in the 
current step by Temephos dissemination. In practice, this requirement seems 
hard to achieve due to several arising technical and spatial heterogeneity 
problems. From Fig. [TJ a suitable choice of constant control pair in the range 
[0,1] x [0,1] still gives insignificant reduction to the number of mosquitoes. 
In an extreme case, killing 100% the number of individuals as in the previous 
step in both Temephos dissemination and fumigation schemes still arises 
endemic situation in the observed area. However, to achieve the condition 
TZ(u) < 1, one has to produce a high value of negative return. For example, 
killing 40% (for indoor larvae and q times indoor eggs) and 180% (for adults) 
as in the previous step will produce 7Z(u) & 0.9. With the same definition, 
killing respectively 60% and 120% leads to 7£(u) < 1 , meaning that there 
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will be a possibility for the mosquito population to be completely eradicated 
for long time. 

Meanwhile from Fig. [2l for whatever the values of (p, /i) are taken in 
the range [0,1] x [1,5], the mosquitoes will never die out. The choice of 
all parameters in the model seems suitable with real situation in endemic 
region, that is, the basic mosquito offspring number with the absence of 
control must be greater than unity. This is another reason why we need 
control intervention. In the figure, /i appears with similar sensitivity as u 2 
in Fig. [H It can be seen that /i is very sensitive, as its slight changes lead 
to the significant differences of the hierarchy of endemicity. In reality, /i 
can be a number that depends on meteorological parameters. It is believed 
that the more environmental condition sustain mosquitoes’ life, the higher 
fjL. The result as in Fig. [2] can be another way to confirm that meteorology 
distribution gives a significant impact to the mosquito abundance in the field. 
The figure also tells us that, based on our model, to where adult mosquitoes 
prevalent to breed is not really a big deal. 

5.2 Time-variant optimal control 

Initially, we consider the trajectory of solution to the equation x = f (x, 0), x 0 ^ 
0, defined as x(t; 0). It can numerically be shown that x(f; 0) is monotonic 
and tends to the coexistence equilibrium Q 4 for very long time. Set ti = 0, 
t '2 — 50, £3 — 100 and 

T>i : = {x : x = f (x, u), t G [0, T], x 0 = x(£j; 0), u G W}, i = 1,2, 3. (22) 

We divide our numerical scheme based on three scenarios: 

(Sc-1) find (x, u) 6 Di x W such that J(u) —> min, 

(Sc-2) find (x, u) G V 2 x U such that J(u) —> min, 

(Sc-3) find (x, u) G x U such that J(u) —> min. 

It is supposed that we want to check the dynamics of mosquito popu¬ 

lation in house-scale including its nearest neighborhood. On average pro¬ 
portion, we let x(fi;0) = [ 8 , 8 , 6 , 6 , 5]. Thus we obtained from our pre¬ 
liminary simulation that x(t 2 ; 0) = [200.41,354.84,73.82,134.85,30.56] and 
x(t 3 ; 0) = [283.12,514.21,91.19,168.22,39.16] respectively. 

One purpose of dividing initial condition for our simulations is to check 
which scenario arising with the cheapest cost, the one which is applied in T 
days after t — t\ (the earliest growth time), t — t 2 (the peak of outbreak) or 
t — tn (the population almost reach the coexistence equilibrium). It may also 
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be the case that the cost increases with respect to the magnitude of initial 
condition. Another purpose is to check which scenario giving the least total 
endpoints ||x(T; u)||i and if the optimal control results in positive value for 
ll x (°;u)||i-||x(T;u)||i, then there is a number n > 1 such that continuation 
of optimal control scheme will completely eradicate mosquito population in 
nT days. Concerning definition of the total cost C(u), we define a weighting 
factor A > 0 such that 



Numerical results from Sc-1 are given as in Figs. [3HH1 




Figure 3: Trajectory of indoor egg Figure 4: Trajectory of outdoor egg 
population. population. 

At first look, it is clear from Figs. [3]|7] that the optimal control makes 
significant reduction to all classes of mosquito population. With the same 
scheme, we can find similar results for both Sc-2 and Sc-3. It turns out that 
by using the data from Table [fl fumigation relatively needs to be enhanced 
rather than Temephos dissemination during application of the control. Figs. [1] 
and [S] tell us so. For further managerial reference, Table [2] represents the per¬ 
formance of optimal control in each scenario. We also highlight the optimal 
control performance when one scheme vanishes. 

From the table, we conclude: 

1 . A complete control intervention using two schemes appears with the 
highest cost for each of three tested scenarios compared to the ones 
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Figure 5: Trajectory of indoor larva Figure 6: Trajectory of outdoor larva 
population. population. 


C*(u), with A = 1 

Sc-1 

Sc-2 

Sc-3 

ui 7 ^ 0, u 2 = 0 


FPPPyffliM 

2.2304 x 10 2 

ui = 0, u 2 ^ 0 

6.6658 x 10 3 


1.0150 x 10 4 

ui ^ 0, u 2 ^ 0 



1.6329 x 10 4 

|| x (0;u)||i, ||x(T;u)||i 

Sc-1 

Sc-2 

Sc-3 

Ul 7 ^ 0, u 2 = 0 

33, 839.57 

794.50, 889.13 

1.09 x 10 3 , 901.57 

ui = 0, u 2 ^ 0 

33, 157.66 

794.50, 154.73 

1.09 x 10 3 , 156.23 

ui ^ 0, u 2 ^ 0 

33, 31.62 

794.50, 35.34 

1.09 x 10 3 , 37.52 


Sc-1 

Sc-2 

Sc-3 

Ui ± 0, u 2 = 0 

610.9374 

946.5915 

11121! • iwsa v ci 

ui = 0, u 2 ^ 0 

196.5743 

225.8798 

235.1042 

ui ^ 0, u 2 ^ 0 

74.7872 

119.5978 

132.6077 


Table 2: The performance of optimal control. 


when one scheme vanishes. It is also noted that the cost slightly in¬ 
creases from Sc-1 to Sc-3, meaning there is a positive correlation be¬ 
tween initial condition taken for simulation and the cost. If we consider 
epidemiology of dengue amongst humans and mosquitoes, this fact can 
help us to decide when we should start to conduct eradication pro¬ 
gramme. It appears as well in the table that there is a huge difference 
between the cost for maintaining only Temephos dissemination and 
only fumigation. Mathematically, this phenomenon arises due to the 
choice of the trade-off coefficient u x . In this case, we assume that the 
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Figure 7: Trajectory of adult popula- Figure 8 : Trajectory of optimal con- 
tion. trol. 

adults need to be accounted for more attention, hence we put a higher 
value for a; X) 5 . However, this does not always mean that taking more 
consideration to reduce adult trajectory requires a higher value of fumi¬ 
gation measure rather than Temephos dissemination measure. In our 
simulation, even if we take w Xj 5 similar to the others, the result is more 
or less similar to the one as in Fig. [8} A reader can play around with 
this coefficient in order to get the result which suits the real situation 
best. 

2. From the table, it is seen that a complete combination of two control 
measures can extremely reduce the total endpoints of mosquito popu¬ 
lation. Maintaining only Temephos dissemination leads us to the worst 
case, i.e. there will be up to 800 total individuals remaining alive after 
T days of treatment. Meanwhile, maintaining only fumigation results 
in moderate total endpoints. If we compute ||x(0;u)||i — ||x(T; u)||i, 
then we would highly recommend to apply both control schemes as an 
integrated programme. 

3. We interpret y Si=ill x *(^> u)||li as ^ ie average number of total pop¬ 
ulation size at each time during observation. The higher its value, 
the higher the number of dengue incidences that may probably occur. 
Again, maintaining only Temephos dissemination leads to the worst 
case. We highlight that the complete control with two measures gives 
the best reduction. Now from the table we know that a preventive act 
of starting the control programme when the number of individuals is 
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small and a simultaneous combination of both control measures will 
bring us tthe best results: least cost, fewest total endpoints and fewest 
number of individuals living during treatment. However, due to practi¬ 
cal limitations, this ideal situation may not be achieved. Knowing the 
density of mosquito population in average house, we can refer to which 
column in the last three columns of Table [2] in order to give the best 
decision for a suitable control management. 


6 Conclusion 

We have an optimal control model for mosquito population dynamics with 
indoor-outdoor life cycle classification. This model initiates the depicting 
situation in most dengue endemic areas where mosquitoes originate from in¬ 
door as well as outdoor. Two schematic control measures are added to the 
model as regulators in reducing the evolution of mosquito population in for¬ 
ward time. The first measure represents a negative return that is earned by 
maintaining Temephos dissemination to some indoor spots where immature 
phases live. Temephos dissemination aims at killing larvae and a small per¬ 
centage of eggs. Meanwhile the second measure is represented by a negative 
return earned from fumigation that targets directly to adults. 

Some brief underlying results from the work in this paper are: 

In constant control case, the proposed model enumerates a primary bi¬ 
ological meaningfulness, i.e. the evolution of mosquito population in for¬ 
ward time results in positive values. With the given estimate of param¬ 
eters and initial condition, it can numerically be seen that the evolution 
of total population remains bounded for all t > 0 by the positive number 
max{||x(0;u)||i, ||Q 4 ||i}. 

If the basic mosquito offspring number is less than unity, then two condi¬ 
tions arise: (i) the zero equilibrium is locally asymptotically stable and (ii) 
the coexistence equilibrium does not exist. 

If the basic mosquito offspring number is greater than unity, then three 
conditions arise: (i) the zero equilibrium is unstable, (ii) the coexistence 
equilibrium exists and (iii) is globally asymptotically stable in nonnegative 
orthant. 

The optimal control results as displayed in Table [2] can be used as a 
reference in decision-making process. A brief conclusion states that: in an 
endemic situation, the best mosquito control impacts are produced if the 
control is started when the number of individuals is as small as possible, also 
with combination of the two simultaneous schemes proposed in this paper. 

Both in constant control and optimal control cases, one has to enhance 
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implementation of fumigation rather than Temephos dissemination. 
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